Josephson charge-phase qubit with radio frequency readout: 
couphng and decoherence 



A. B. Zoriifl 

Phystkahsch-Technische Bundesanstalt, Bundesallee 100, 38116 Braunschweig, Germany 

(Dated: 9 December, 2003) 

The Josephson charge-phase qubit based on a superconducting single charge transistor inserted 
in a low-inductance superconducting loop is considered. The loop is inductively coupled to a radio- 
frequency driven tank circuit enabling the readout of the qubit states by measuring the effective 
Josephson inductance of the transistor. The effect of qubit dephasing and relaxation due to electric 
and magnetic control lines as well as the measuring system is evaluated. Recommendations for the 
qubit operation in the magic points with minimum decoherence are given. 
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INTRODUCTION 

The superconducting quantum bit (qubit) circuits 
comprising mesoscopic Josephson tunnel junctions have 
recently demonstrated remarkable quantum coherence 
properties and now are considered to be promising el- 
ements for a scalable quantum computer The read- 
out of macroscopic quantum states of a single qubit or 
a system of coupled qubits with minimum decoherence 
being caused by the detector, remains, however, one of 
the most important engineering issues in this field. 

The Josephson qubits are commonly subdivided into 
flux, phase, charge and charge-phase qubits. The design 
of charge and charge-phase qubits is based on a Cooper 
pair box in which a small superconducting island with 
significant Coulomb energy is charged through a small 
Josephson junction (charge qubit) or a miniature double- 
junction SQUID (charge-phase). The distinct quantum 
states of the box generated by signals applied to a gate 
are associated with different observable charges on the 
island. This makes it possible to read out the qubit 
state by discriminating the island charge. Probing this 
charge can be done either by single quasiparticle tunnel- 
ing across a small auxiliary tunnel junction attached to 
the island or by a capacitively coupled electrometer 
In the charge-phase qubits, the quantum states of 
the box involve the phase coordinate of the SQUID loop, 
so discriminating these states can also be done by mea- 
suring the persistent current circulating in the loop at 
an appropriate dc flux bias. Such a measurement was 
performed in the experiment of the Saclay group 
their setup, nicknamed " Quantronium" , the circulating 
current passed through a larger auxiliary (third) junction 
was read out by measuring the switching current of this 
junction. 

The persistent current is not, however, the only phase- 
dependent quantity characterizing the quantum state of 
the charge-phase qubit. Another useful quantity is the 
Josephson inductance of the double junction, which can 
be probed by small rf oscillations induced in the qubit. 
Recently, we proposed a transistor configuration of the 



Cooper pair box (see Fig.^, in which the macroscopic 
superconducting loop closing the transistor terminals was 
inductively coupled to a radio- frequency tank circuit . 
Similar to the rf-SQUID-based method of measurement 
of the Josephson junction impedance 0, this setup makes 
it possible to measure the rf impedance (more specifically, 
the Josephson inductance) of the system of two small 
tunnel junctions connected in series and in doing so, to 
probe the macroscopic states of the qubit. 

On the one hand, the advantage of this method con- 
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FIG. 1: (a) The electric circuit diagram of the charge-flux 
qubit inductively coupled to a tank circuit by mutual in- 
ductance M . The macroscopic superconducting loop of in- 
ductance L is interrupted by two small Josephson tunnel 
junctions positioned close to each other and forming a sin- 
gle charge transistor; the capacitively coupled gate polarizes 
the island of this transistor. The qubit is controlled by charge 
Qo generated by the gate and flux $m induced by coil Lm- 
The tank circuit which is either of a parallel (b) or a serial (c) 
type is driven by a harmonic signal (J^f or V^t, respectively) of 
frequency tUrf ~ ujq, the resonant frequency of the uncoupled 
tank circuit. 
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sists in an effective decoupling between the qubit and a 
measurement device, tliat reduces the decoherence of the 
qubit. Moreover, the loop design of the qubit has a po- 
tential to perform a data readout in a non-destructive 
way 0. On the other hand, due to the selective charac- 
teristic of the tank, the bandwidth of this setup is rather 
narrow, so the optimum relation between the relaxation 
time of the qubit and the time of measurement becomes 
an issue. Furthermore, the driving rf signal may induce 
appreciable frequency modulation and dephasing of the 
qubit during its evolution (performance of the quantum 
operations). Switching the oscillations off and on is, how- 
ever, possible only on a relatively long-time scale of a 
transient process in the tank. 

In this paper we address the problem of decoherence 
induced in the charge-flux qubit by the classical reso- 
nance tank circuit. Besides, we propose a strategy of 
measurement and optimize the regime of qubit operation 
for typical parameters of the circuit. 



BACKGROUND 

The small tunnel junctions of the charge-flux qubit 
are characterized by self-capacitances Ci and C2 and 
the Josephson couphng strengths Eji and Ej2- These 
junctions with a small central island in-between and a 
capacitively coupled gate therefore form a single charge 
transistor connected in our network as the Cooper pair 
box (see Fig.^). Critical currents of the junctions are 
equal to /ci,c2 = ^Ejij2, where $0 ~ h/2e is the flux 
quantum, and their mean value Ico = ^{Ici + Ic2)- The 
design enables magnetic control of the Josephson cou- 
pling in the box in a dc SQUID manner. The system 
therefore has two parameters, the total Josephson phase 
across the two junctions (j) = + (p2 = 27r$/$o con- 
trolled by flux $ threading the loop and the gate charge 
Qo set by the gate voltage Vg. The geometrical induc- 
tance of loop L is assumed to be much smaller than the 
Josephson inductance of the junctions Ljq = $o/(27r/co). 



/3l = L/Ljo < 1. 



(1) 



Neglecting the magnetic energy term associated with 
the current through the small inductance L, the Hamil- 
tonian of the autonomous qubit circuit is expressed as 



\Eji-Ej2\ < Ej{cj)) < E.n+Ej2 = 2Ejo = ^I^o, while 
a phase variable x = + 7(0)- The angle 7 is given by 
the expression 



tan7 = (ji - j2)tan(0/2). 



(4) 



where dimensionless Josephson energies ji,2 — 
£'Ji,J2/(2i?jo) with ji + j2 = 1. The difference 
phase (p = ■^{(fii — f2) is a variable conjugate to the 
island charge 2en = —2ei-^ and n is the operator 
of the number of excess Cooper pairs on the island. 
This charge enters the charging energy (first) term 
in Eq. (O, in which C is the total capacitance of the 
island, C ~ Ci + C2 + Cg ~ Ci + C2 and the gate 
capacitance Cg ^ Ci.2- The characteristic charging 
energy, Ec = e^/2C is assumed to be of the order of the 
Josephson coupling energies Eji Ej2 3> ksT. 

The Schrodinger equation correspondingto the Hamil- 
tonian Eq. ^ is the Mathieu's equation l9|. The eigen 
energies form the Bloch bands and eigenfunctions \n, q) 
are the Bloch wave functions of a particle in the peri- 
odic (Josephson) potential with " quasimomentum" (here 
rather quasicharge) q. Its value is the charge which the 
gate source provides to the island, i.e., q = Qo = CgVg. 
Each of such eigen functions can be presented as a coher- 
ent superposition of the plain waves. 
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where m = 0, ±1, ±2, ... is the number of excess Cooper 
pairs on the island |lClllll| . The weights of these coherent 
contributions |Cm^p depend on q, the band index n and 
the characteristic ratio 



X^Eji^)/E, 



(6) 



The two lowest energy levels En{q,4>), i-e., n = and 1 
(see their dependencies on q and (f> in Fig.lJJ, form the 
basis {|0), |1)} suitable for qubit operation. In this basis, 
the Hamiltonian Q is diagonal. 
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(7) 



where ai with i = x,y,z is Pauli spin operator. The 
general state of qubit is 



(8) 



Ho^^^^^^^~EA^)cosx. (2) 

The second term in Eq. (0) originates from the total 
Josephson energy equal to —Eji cos (pi — Ej2 cos ip2 ■ The 
effective Josephson coupling strength is 



Ej{4>) = {Ej^ + E% + 2E.nEj2 cose 



,1/2 



(3) 



with 



|6p = 1. It is remarkable that the level spac- 



ing e{q,<j)) = hn = Ei{q,(jji) — Eo{q,(j)) and, therefore, 
the transition frequency f2 are efficiently controlled by 
two knobs, i.e., by varying the parameters q and (or, 
equivalently, Qq and $) (121] . 

The idea underlying the measurement of this charge- 
flux qubit is based on inducing radio frequeiicv oscilla- 
tions in the tank circuit of frequency Lu^f <IC [g . Due to 
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inductive coupling M, these oscillations cause oscillations 
of corresponding flux <i>T (see Fig. la) and, therefore, of 
total phase. 



27r 



($T + $,„) =0aSin(wrft + 6l)+(/.o. (9) 



If the rf drive signal is sufficiently weak, the amplitude (pa 
of these oscillations is relatively small, pa <C tt. In this 
linear regime the reverse Josephson inductance is equal 
to 



(10) 



i.e. it is determined by the local curvature of the energy 
surface £"„. For example, for Ejq = 2Ec (see Fig. 2) at 
q w 0, we have the following estimates within the zero 
and first bands: 



and 



(0,0,0) « ALj^ cos p, 



Lj\l,0,(j)) « O.liJo cos 



(11) 



(12) 



respectively. In the vicinity of the level avoided crossing 
point, q = e and (p — tt {in Fig. 2 marked as D-D'), the 
reverse inductances may increase significantly, 

L7\n,0,7r) ^ \ . . , L-j},, n = and 1, (13) 

due to small asymmetry of the transistor parameters, 
\ji — j2\ 1. For example, in the case presented in 
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FIG. 2: Shape of energy bands £"0 and Ei in the charge- 
flux qubit calculated for the mean Josephson coupling Ejo = 
^{Eji + Ej2) = 2Ec and the Josephson coupling asymmetry 
parameter \ji - j2 \ = {Eji ~ Ej2)l{Ej\ + Ej2) = 0.1. Black 
(hoUow) circles on the zero (excited) band surface mark the 
locations of magic points A {A'\ B {B'), and C (C') and the 
energy level avoided crossing point D (D'). 



Fig. 2, |ji - J2I = 0.1 and L^^ = =F2.5Ljp for the zero 
and first band, respectively. In points C and C", the ab- 
solute values 1^7^! are smaller but the signs for n = 
and 1 are still different. 

Coupling to the qubit causes a shift of the resonance 
frequency t^o = {LtCt)~^^^ of the tank circuit, i.e., 
^o("-) — '^o + S^oin), where 

SLOo{n)^-h^(3L , ,, u;o. (14) 



' Lj{n,q,(l)) 

Here k is the dimcnsionless coupling coefficient, 

M 

< 1. 



(15) 



The resonance frequency shift duJo{n) carrying informa- 
tion about the qubit state \n) is found from the ampli- 
tude or/and phase of forced oscillations in the tank. For 
achieving sufficient resolution in such measurements, the 
quality factor of the tank circuit Q should be about, or 
larger than, the ratio ujo/\Sujo{Q) — '^^o(l)|- 



INHERENT AND EXTERNAL SOURCES OF 
DECOHERENCE 

In our consideration we had neglected the quasiparticle 
tunneling which inevitably causes dissipation of energy. 
Even rare tunneling of individual quasiparticles across 
the tunnel junctions, i.e., on and from the island, can 
decohere the qubit and completely destroy the readout 
regime described above. These processes lead to sudden 
change of the operation point, q —> q ± e and, possibly, 
the energy band index, i.e., cause relaxation 1 — > 0. 

The processes of single quasiparticle tunneling across 
a small Josephson junction had been studied by Averin 
and Likharev in Refs. 0, ■ They had generalized the 
orthodox theory of single electron tunneling to the case 
of finite Josephson coupling, Ej ^ 0, taking into account 
the dynamics of essential phase factors exp(±ix/2) in 
the electron tunneling terms added to the Hamiltonian 
of type Eq. (|2Jl . These factors are the operators of single- 
electron transfer and their nonzero matrix elements in 
our basis are 

^nn' = {n,q\exp{±ix/2)\q±e,n'). (16) 

The rates of transitions \q, n) \q ± e, n') are given by 

-r± _ u± 1 2 -^qp(^K7i'/s) 



1 — exp — 



(17) 



In our case Iqp(U) is the quasiparticle current-voltage 
dependence of the network of two tunnel junctions of 
the qubit connected in parallel. Due to 2e-periodicity of 
the energy surfaces, the corresponding energy gains are 
identical. 



En{q,P) - En>{q±e, 



(18) 
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FIG. 3: Off-diagonal matrix elements of the single quasiparti- 
cle transfer operators exp(±jx/2) computed for different val- 
ues of equivalent Josephson coupling set by flux <I> = $0/2 
(solid lines), $o/4 (dashed lines) and (dotted lines). The 
qubit parameters are the same as in Fig. 2. 



and their value depends on the operation point {Qo: 
(see Fig. 2). 

The relation between this energy and the supercon- 
ductor energy gap Age is important for having the quasi- 
particle transitions infrequent or even get rid of them. 
First, if voltage U = e^„//e < 2Asc/e, the quasipar- 
ticle current Iqp{U) entering Eq. ()17|l is exponentially 
small, i.e., ~ /co exp(— Asc/fcs?") 0|- At larger volt- 
ages, U > 2Asc/e, current Igp is enormously large, 
> 2/co. Therefore, in order to prevent intensive tun- 
neling of quasiparticles the gain of energy, e^^, , should 
never exceed 2 Age. Secondly, if this gain is smaller than 
Age, then infrequent quasiparticle tunneling can, in prin- 
ciple, be blocked by the gap energy associated with one 
unpaired electron in the superconducting island (the so- 
called even-odd parity effect) p^ . 

Suppression of quasiparticle transitions within the zero 
energy band in superconducting Al single charge transis- 
tors and Cooper pair boxes was extensively investigated 
experimentally. Dep endin g on experimental skill and 
luck (see, e.g., |l3)llMllMl23) inspected devices often 
exhibited pure Cooper pair behavior when their charg- 
ing energy Ec was not larger than ~ 100 /zeV « O.SAai, 
where Aai is superconductor energy gap of aluminium. 
Since the energy gain for transitions in the Cooper pair 
boxes and low- voltage-biased transistors, Cqq, is less than 
Ec for any Ej, the condition Ec < Age can in a "good" 
qubit sample ensure suppression of quasiparticle tunnel- 
ing in the ground state. 

For quasiparticle transitions from the excited state this 
condition is clearly insufficient. For example, for small 
Ej (corresponding to flux value <I> = <I>o/2, Eq. ((SJ) the 
energy gain values are between about Ec (for the process 
D' A, see Fig. 2) and 4Ec (for the processes A' D' 
and A' —^ D). At larger Ej both the minimum and max- 
imum energy gain values are even larger. For example, 
for Ej = AEc (i.e., $ = 0), transitions C B and 
B' ^ C correspond to energy w AEc and ~ 5Ec, respec- 
tively. Since the first multiplier factor in the expression 



for resulting relaxation rate, 

[T<^^P^ = rjo + rro^ (k^oP + kroH ^^'^^f'^'^ (19) 

is nonzero for any Qq and <& (see the plots of the two 
items in Fig.O, only the condition Ec < Asc/5 can en- 
sure suppression of these transitions in arbitrary oper- 
ation point of our qubit. Possibly insufficiently small 
value of Ec was the reason of very short relaxation time 
(tens of ns) in the recent experiment with a charge qubit 
by Duty et al. [i^. Their Al Cooper pair box had 
Ec ~ 0.8 Age and Ej « 0.4i?c, so the energy gain in the 
chosen operation point (Qq — 0.4e) was too large, i.e., 
about 2.2Ec ~ 1.8 Age > Age (although in the ground 
state this sample nicely showed the pure Cooper pair 
characteristic) . 

Moreover, there are several sources of decoherence due 
to coupling of the qubit to the environmental degrees 
of freedom. For evaluating the effect of these sources 
on the qubit the coupling Hamiltonian term, -ffcoupi = 
njf'^ + Hc"\ is included in the total Hamiltonian of the 
system, 

-ff = -ffo + -ffcoupl + -H^bath , (20) 

where i/bath is a bath operator. Hc"^^ and Hc"^^ are the 
electric control line term and the magnetic coupling term 
respectively. The latter is associated with both the flux 
control line and the tank circuit. Fluctuations originating 
from the sources of gate- and flux-control lines may, in 
principle, lead to significant decoherence of the qubit. As 
was shown in Ref. |2ll | and demonstrated in experiments 
0, 0, 0, these effects can, however, be minimized by 
choosing appropriate (minimum) coupling. On the other 
hand, decoherence caused by the tank-circuit-based read- 
out system requires special analysis, because weakening 
of this coupling leads the input signal being reduced. Be- 
low we will start with sources of decoherence associated 
with the control lines and then will analyze the effect of 
the tank circuit and amplifier. 

COUPLING TO CHARGE CONTROL LINE 

The coupling of the charge-phase qubit to the elec- 
trical control line is actually similar to that of the gate 
coupling in the ordinary Cooper pair box 0]. However, 
here we assume that Josephson coupling parameter A is 
not necessarily small, as is usually assumed in the analy- 
sis of the charge qubits. This generalization of the model 
is essential because the external flux $„i changes the ef- 
fective Josephson energy Eq. (j^J of the qubit over a wide 
range. The assumption of not small A implies that the 
eigenstates of our system Eq. ^ are generally composed 
of several (not only two) plane-wave states. 
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The coupling term can be presented as 

H^/^ = -2enSVe, 



(21) 



where 5Ve is the operator of vohage fluctuations on the 
island in the absence of Josephson coupling. The charge 
operator is equal to 2en — Qq~ CV, so the essential part 
of the coupling Haniitonian is 



The voltage operator is equal to 

,> $0 . $0 / . dj 



$0 

2tt 



(22) 



(23) 



Here we assumed slow variation of the total phase 0, 
Eq. iPJ. The operator of voltage V is similar to the op- 
erator of velocity of an electron in the periodic electric 
potential of the crystal lattice , and its interband ma- 
trix elements are 



-dn,n' +1 Xnn'{l - On,n'), (24) 



dq 



2e 



where 5n.n' is the Kroneker delta and Xnn' are the matrix 
elements of the phase operator x ^3 ■ 

Finally, the coupling Hamiltonian, Eq. H22|l . takes the 
form 

H'f^ = {ax sin rje + <Jz cos r\e)Xe-, (25) 
where operator 

X,^C\\V\\We (26) 
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FIG. 4: The terms composed of diagonal (a) and ofF-diagonal 
(b) matrix elements of operator V and entering Eqs. <\27^ and 
H28|l are presented for different values of flux <l?e for the given 
qubit parameters (see caption of Fig. 2) . 



with 



and 



\\V\\ = -V(^4i-K)o)2+4|l-oiP 



tan ?7e 



2\Voi\ 



(27) 



(28) 



(The plots of the terms entering Eqs. (|?7|) and (HHJ ob- 
tained by numerical calculations are presented in Fig.0]) 
Thus, Xe =^2a ^aXa Can bc considered as an operator 
of the bath l23 with the Hamiltonian 



E 



2ma 



(29) 



and the spectral density S^^' i"^) ^ C^WW^ s'^\^) = 
Je{uj)Q{uj,T)/iLj. Here the oscillator energy function is 



and 



(30) 
(31) 

Assuming smallness of Cg, the spectral density of 
the fluctuations of SVe is given by 

4'^(^) = ^(^) ReZt{uj)Q{Lo,n (32) 

where Zt = (iujC'g+Zg^)^^ is determined by parallel con- 
nection of the qubit capacitance CgC/ {Cg + C) ^ Cg and 
the gate line impedance ZgiiS) ~ i?ioo = 100 17. There- 
fore, for frequencies up to LOg = (i?iooCg)~"'^ ^ e/^j i-e-j 
at all characteristic frequencies of the system, ReZt = 
i?ioo • This is the case of linear damping in the Caldeira- 
Leggett model. 



with the dimensionless factor 

/ r , _ 1 1 1/ 1 1 \ 



c,\\v\\ 



Ri 



< 



R 



100 



Rc 



(33) 



(34) 



where Rq = h/Ae^ ~ 6.45 kfi, the resistance quantum. 
The estimate similar to the last expression in Eq. I|34(l . 
was given earlier in Ref. [Ui for small A. 

Relaxation and dephasing caused by the charge control 
line can, therefore, be described by the spin-boson model 
with linear damping The corresponding rates are 

given by the expressions 



TTCke sin rje ^ coth 



2kBT' 



(35) 
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and 



One can see that in accordance with the conclusions of 
Refs. 0, 113, reducing the coupHng coefficient ag by a 
small factor (Cg/C)^ ^ 1 can significantly depress the 
decoherence rates. 



COUPLING TO FLUX CONTROL LINE 

The inductive coupling of the qubit loop with the con- 
trol and readout circuits is described by the Hamiltonian 



where Is is the operator of current circulating in the qubit 
loop, d^rn = MmSIrm the bath operator (proportional 
to fluctuations of current Sim in the control inductance 
Lm); S^T = MSI is the operator of flux associated with 
current fluctuations in the tank circuit. 

To specify the coupling, we present the operator Is in 
the eigen basis (jSJ, i.e., we find the matrix elements 



{n\is\n'), n,n' = 0, 1. 
In the general case. Is is given by the expression 

is = KlQl + K2Q2, 



(38) 



(39) 
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FIG. 5: Terms composed of diagonal (a) and off-diagonal (b) 
matrix elements of operators cos x and sin x, respectively, cal- 
culated for different values of flux $e for the given qubit pa- 
rameters (see caption of Fig. 2). 



where the dimensionless factors K1.2 
identity K1 + K2 = 1- Quantities Qi,2 - 



C2,i/C with the 
— 2e?75-^ — are the 



(36) charges on the first and second junction, respectively, and 
their time derivatives are the Josephson supercurrents. 



Qi,2 = t[Qi,2,H()] = /ci,c2sin(^i 2. (40) 

Using identity 1^91,2 = f±V' = 2^^"F7 and Eq. the 
circulating current can be presented as 



is = /i (0) cos X + /2 (0) sin X. 

The amplitudes of these two components are 

2nE,nEj2 . . 
Ii = — „ , smc 



(37) and 



(jl - J2)('«ljl +K2J2): — 



(41) 



(42) 



4tt EjiEj2 2 
+ ^'^^"'^^^^^7(^'°^ 2' 



(43) 



respectively. 

Since the Hamiltonian, Eq. Q, is an even function of 
X, the operators cos x and sin x entering Eq. (|41|l are di- 
agonal and off-diagonal, respectively. Amplitude Ii is 
merely the classical Josephson current across two large- 
capacitance junctions as a function of the overall phase 
difference 4>, while the diagonal term cos x describes the 
suppression of this current due to the charging effect 
{Ec 0) (see, for example, Ref. 0|). The second, off- 
diagonal, term in Eq. (|41|l is due to asymmetry of the 
transistor; it gives rise to the interband transitions <-*■ 1. 
Denoting the nonzero values of the corresponding matrix 
elements as 



Coo 



|cosx|0), cii = (1| cosx|l), 



(44) 



and 



SOI = |(0|sinx|l)|, (45) 
we obtain the coupling Hamiltonian Eq. H37|) in the form 

Hi'"^ = {crv sin 77™ -f cos 77™ ) (X„ + Xt ) , (46) 
where 

Xm,T^-\\I\\S<P,n,T, (47) 



= V[(cii-coo)/i]2 + [2soi/2]2, 



tan?7j, 



2S01J2 
(cii - coo)/i ' 



(48) 



(49) 
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(see the plots of terms entering these expressions in 
Fig.El). 

First we omit in Eq. H4()|) the term Xx associated with 
fluctuations of the tank circuit and focus on the effect of 
fluctuations in the flux control line (5$m = MmSIm- As- 
suming real impedance of the flux control line, Z,„ ^ 
i?ioo, we obtain the spectral density of the operator 



oc 5Im in the form S^;^\io) = MlWIfSr'ii^) 
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Jm(w)0(w, T)/a;. At frequencies below uj„ 
the function J^, is linear, 



Rim/ Ln 



(50) 



with the dimensionless coupling factor 

/ "/ /l/7.__ II / II \ ^ 



UJ > UJj, 



Rq 



R 



100 



(51) 



the effective damping 



At higher frequencies, 
decays as (cj,„/a>)^. 

In fact, Eq. (|51|l describes the effect of coupling to the 
control flux in the general case. An estimate of the cou- 
pling factor based on the evaluation 



l/ll 



dEj 



9$ 



(52) 



which is valid for symmetric transistor {I2 = 0) with 
small Ej, was made by Makhlin et al. Small mutual 
inductance Mm ^21,. 25] leads to small a™ and, therefore, 
causes significant depression of the corresponding relax- 
ation rate. 



sin^ ?7„ 



rt coth 



nn 

2kBT' 



and dcphasing rate. 



-1 ^ 2 ^ksT 

+ na,n cos r]m — - — . 



(53) 



(54) 



DECOHERENCE DUE TO READOUT SYSTEM 

In contrast to the control lines, coupling to a read- 
out device (in our case the tank circuit with amplifier) 
cannot be made arbitrarily small in order to reduce the 
decoherence. This coupling should ensure sufficiently 
strong signals at the amplifier input in order to perform a 
measurement with a reasonable signal-to-noise ratio on a 
time scale shorter than that determined by other factors, 
namely Tr'^K Moreover, without an efficient switch (see 
possible design of such a switch in, e.g., Ref. 27]), such 
a coupling may cause significant dephasing of the qubit 
during quantum gate manipulation. 

Inductive qubit coupling to the tank circuit is de- 
scribed by the Hamiltonian Eq. I|46|l . The spectral den- 
sity of fluctuations of the corresponding variable Xt oc 
6^T = MSI is expressed as 



TT 



,H.)2<-i:i. (57, 



where Sj {uj) is the spectral density of the noise current 
SI across inductance Lt- Since the cold (superconduct- 
ing) tank circuit itself presumably has very low losses, a 
backaction noise SI of the amplifier is dominating. It is 
associated with input real impedance, modeled by Rp or 
Rs for parallel and serial configurations, respectively (see 
Fig. 1). T* is the effective temperature associated with 
this impedance. 

Spectral density Sj and function Jt{oj) can be found 
from a network consideration. Neglecting small detuning 
SujQ <^ luq, in the case of the parallel network (Fig. lb), 
the spectral function Jt is given by the expression 



^"f"~(t^2_^2)2_,.^2^2g-2' 



(58) 



with 



So far we have considered the effects of decoherence 
due to the charge and flux control lines as two indepen- 
dent effects. Actually, one has to describe them together 
using a multibath model [IJ. If either of these deco- 
herence effects is small, i.e. the so-called Hamiltonian- 
dominated regime is realized, the total rates due to con- 
tributions of the two control lines are given by 



(c)l 



r(™)l 



(55) 
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cO 



(59) 



where the quality factor Q — loqCtRp — Rp/^qLt- 
the case of the serial network shown in Fig. Ic, 



In 



with 



-asIiLo 



2 , ,2\2 



,2,. ,2/ 



(60) 



lrW]-i ^ [2rW]-i ^ [^(e)]-i ^ 



(56) 



In our model we shall assume that such a regime is re- 
alized and, moreover, the resulting rates, Eqs. H55(l and 
()56|l . can be made negligibly small. Below we will fo- 
cus on the effect of the readout circuit, whose coupling 
strength has to be optimized. 



as 



f2Am 



RqRs 
(^^ 



(61) 



and Q = {ujqCtRs) ^ = ujqLt/Rs- 

In contrast to the linear spectral functions for the 
control lines, Eqs. and (|5n|l . the functions given by 
Eqs. (|58|l and (|60|l describe a structured bath, viz. they 
both are of a Lorentzian (resonance) shape. A similar 
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situation emerges, for example, in the case of the flux 
qubit with readout using a C-shunted dc SQUID [isf . 
The theoretical analysis of the spin-boson model with a 
structured bath was made by Kleff et al. on the ba- 
sis of the flow equations. If coupling is, as in our case, 
weak, only the high frequency (w ^ 57) and low frequency 
{uj —^ 0) behaviors of J (a;) account for relaxation and de- 
phasing, respectively [j^Llifll Isilj l. 

Since frequency Q is typically about tens of GHz and 
the distance between qubit and amplifier presumably ex- 
ceeds the wavelength, the effective real admittance of the 
parallel circuit at these frequencies is equal to i?j~gQ and 
the impedance of the serial circuit « i?ioo- Therefore, the 
relaxation rates increase by factors gp = Rp/RiQo ^ 1 
and gs = Rioo/Rs ^ 1, respectively. 

For the parallel tank circuit the relaxation and dephas- 
ing rates (presumably, <C ujq) are equal to 



Q, coth 



nn 



and 



[r(P)]-i = [2r(P)]-i+^apCOs2r?„ 



2kBT* 



(62) 



(63) 



respectively. The relaxation rate is dramatically sup- 
pressed due to the small frequency ratio, {luq/H) <C 1. 
For the serial configuration, the corresponding rates are 



V n 



2^ ^ nn 

il coth 



2kBT* ' 



(64) 



(65) 



The dephasing rate is determined by the rate of re- 
laxation, because at low frequency, w ^ wq, function 
cx UJ j2J|. Due to weaker decay of the serial 
circuit impedance at high frequencies, lo ^ loq, the re- 
laxation rate is, however, substantially higher than that 
in the case of the parallel circuit configuration. We shall 
therefore focus our further consideration only on the par- 
allel tank circuit as the more favorable (allowing longer 
measuring time). 



MAGIC POINTS AND SOME ESTIMATIONS 

The analysis of the coupling between the qubit and 
the tank circuit, Eas. and Fig. (5), shows that 

its strength Xt oc \\I\\ and mixing angle 77™ can be signifi- 
cantly varied by choosing an appropriate operation point. 
For example, as can be seen from Eq. (|42l) . the diagonal 
component of Xt (oc Ii) which essentially cause pure de- 
phasing of the qubit, is zero, i.e., mixing angle 77™ — 7r/2, 
at phase values cj) = and tt. The derivatives and. 



therefore, the circulating supercurrent are zero. More- 
over, as illustrated in Fig. 5b, if the gate charge Qo ~ 
(i.e., derivatives ^gq^ = 0), then |soi| and, hence, Xt 
are minimum. In particular, at the bias flux — 
or, equivalently, 4> = tt (this point is marked in Fig. 2 as 
A), Ej{<f>) = \Eji — Ej2\ <C Ec, so one can use the ex- 
plicit expressions for the wave functions, Eqs. (All) and 
(A12) of Ref. [13, and obtain 



\sqi\ 



16%/2 Ec 
Then the value of ||/|| gHl) is 



1 Ej{^) _ \j1-j2\Ejo 
8%/2 Ec ' 



1 7-11 01 1 7- ~ I ji - j2 I Ejo 

|/|U = 2|.oi|/2«^^-^/cO, 



(66) 



(67) 



where we have taken into account that Tti ^ K2 ~ 0.5 
and the second term in Eq. H43|) . cx (ki — K2), van- 
ishes due to the relation cos((/)/2) — 0. In the point 
Qo = 0, (/) = (marked as B in Fig. 2), the Josephson 
energy Ej{(f)) — 2Ejo and 2|soi| is approximately equal 
to {l/8V2)Ejo/Ec, so 



l/ll 



\jl - .h + Kl- K2\ E, 



JO 



IcO, 



8V2 Ec 
while for Qo — e (point C in Fig. 2), |soi| ~ 0.5 and 



(68) 



l/ll 



(69) 



It is remarkable that the effect of asymmetry in critical 
currents and capacitances of the junctions can, in princi- 
ple, cancel if (ji — ^2) — ~{ki — K2). In practice, however, 
the signs of {ji — J2) and (ki — K2) normally are similar 
because the critical current and capacitance are both pro- 
portional to the junction area and such cancelling does 
not occur. 

Comparing Eqs. l(H7|) . and one can see that 
on the assumption of small asymmetry of the transistor, 
ji ~ j2 ~ Kl « K2 ~ 0.5, the coupling strength Up 
in points A{Qo = 0, = tt), B {Qq = 0, </> = 0) and 

(Qo = e, = 0) is rather small, but it is significant in 
the point, D (Qo = e, </> = tt), where parameter |soi| « 
0.5 and 



(70) 



To illustrate this behavior, the coupling strength evalu- 
ated for typical parameters of the system is presented in 
Tablet. 

From the point of view of operation with a minimum 
dephasing rate, the "magic" points A, B and C, in which 
supercurrent Ii — (see Eq. (03)), are clearly preferable 
to other points in the Qq-^ plane. Therefore, manipula- 
tion of the qubit can, in principle, be performed in any of 
these points. The estimated values of the corresponding 



quality factor for quantum manipulation, Qu, = flr^ 



(p) 
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TABLE I: Evaluated qubit parameters derived on the assump- 
tion Ejo = 2E^ = 80/ieV (i.e., ho ~ 40 nA and SS^ = Aai ~ 
200 /ieV, tiie energy gap of Al) and ji — j2 — k\ — K2 = 0.1. 
The tank circuit quahty factor Q = 100, frequency uoq — 
2tv X 100 MHz, (Lt/CtY^'^ = IQOn, k^QPL = 20 and tem- 
perature T* = 1 K 3> r ~ 20 mK. As long as the dephasing 
rate in the magic points is nominally zero, a 0.1% inaccuracy 
of the adjustment of the values (j> = n and was assumed. 



Operation point: 


A-A' 


B-B' 


C-C 


D-D' 


Frequency Q./2-k [GHz] 


39 


50 


36 


4 


Coupling strength Ofp 


210~^ 


10^2 


410-^ 


1 


Qubit quality factor 


310** 


210^ 


lO"! 


< 30 


Relaxation time ri-^^ [s] 


810"^ 


10-1 


610"^ 


10-^ 



given in Table 1, are sufficiently high. For example, in the 
case of preparation of the qubit in point A, the manipula- 
tion can be performed by means of a dc pulse applied to 
the transistor gate dllll^. This pulse (with short rise 
and fall times) can rapidly switch the qubit, for example, 
to point D and back to A causing its evolution (although 
with significant dephasing) during the pulse span. Our 
qubit prepared in the ground state in point A or i? or C 
can be (preferably) manipulated by a pulse of microwave 
frequency, ~ il, applied to the gate. For example, the 
Quantronium qubit in the experiment by Vion et al. j5j] 
was manipulated by microwave pulses in point C. 

For reading out the final state, the qubit dephasing 
is of minor importance, while the requirement of a suf- 
ficiently long relaxation time is decisive. Moreover, the 
relaxation rate may somewhat increase due to oscillations 
in the tank induced by a drive pulse (see Fig.|3, which 
leads to the development of oscillations around a magic 
point along (p axis, Eq. (|5J|. If the frequency of these os- 
cillations is sufficiently low, Wrf ^ they result only 
in a slow modulation of transition frequency r2. Increase 
in amplitude of steady oscillations up to (/>a ~ 7r/2 (de- 
termined by the amplitude of the drive pulse and detun- 
ing) yields a large output signal and still ensures the re- 
quired resolution in the measurement provided the prod- 
uct k^Q^L > 1 is sufficiently large. (At larger amplitudes 
00, the circuit operates in a non- linear regime probing 
an averaged reverse inductance of the qubit whose value, 
as well as the produced frequency shift Sujq, is smaller 
[sif.l As points A and B lie on the axis Qo = and 
are both characterized by a sufficiently long relaxation 
time, reading-out of the qubit state with the rf oscilla- 
tion span ±7r/2 is preferable in either point. In the case 
of operation point C, the limited amplitude of the oscil- 
lations does not reduce much the relaxation time either. 
Significant reduction of the relaxation time occurs in the 
vicinity of point D. Due to this property which is due to 
the dependence of the transversal coupling strength on 
(/), Eqs. the measurement of the Quantronium 

state using a switching current technique was possible in 
the middle of segment CD (see Fig. 2), where the maxi- 



mum values of the circulating current in the excited and 
ground states were of different sign Q . 

In the vicinity of energy level avoided crossing point 
D, in which the gap between the zero and first excited 
states is minimum, hVl = 2\ji — j2\Ejo, oscillations of 
(j) may cause Landau-Zener transitions |0) <-* |1) [s^. 
Probability of a such transition per single sweep 



Plz = exp 



-27r 



(71) 



can be appreciable in a sufficiently symmetric transis- 
tor and/or at high driving frequency to^f, i.e., when 
bi — j2\ ^ ifi^ii/ EjqY^'^ . These transitions lead to un- 
wanted mixing of the qubit states [33l |. In the vicinity 
of point A', where the gap between the first and sec- 
ond (not shown in Fig. 2) energy bands is smaller |^, 
MI12 = (ji — j2)'^Ejf^/2Ec, the Landau-Zener transitions 
|1) ^ |2) are more probable. Fortunately, the second en- 
ergy band has the positive curvature, ^ ^^g^'a'^"^'' > 0; 
the mixing of these states might even improve the distin- 
guishability of signals from the ground and excited states. 
More rigorous analysis of this effect on operation of the 
qubit in point A is, however, needed. 

Finally, let us evaluate the time of measurement re- 
quired for the resolution of the states n — and 71 = 1 in 
the most favorable magic points A and B. As schemati- 
cally shown in Fig. 6, an rf drive pulse is applied to the 
tank circuit just after manipulation of the qubit {t — 0) 
and induces growing oscillations in the tank. The am- 
plitude of the oscillations of voltage V approaches a 
steady value Aq (Ai) for n = Q {n = 1) after the time 
trisc ~ 2TrQ/ujQ. Assuming a corresponding amplitude of 
oscillations of phase (j)a = 7''/2, we obtain for the param- 
eters of Table 1 the amplitudes 



^0 = (I), 



$0 loLt _ f n<i>oUjRpIc 
''2tt M ~ \ Sk^QPL 



1/2 



3 /iV (72) 



and Al w 1 ^V. 

Assuming that the equivalent noise of a semiconductor- 
based amplifier referred to the input is of the order of 
Johnson voltage noise across Rp ^ 10 kO at ambient 

temperature T* ^ 2K, i.e., Sy^ ~ lnV/\/Hz, we can 
express the signal-to-noise ratio as 



SNR: 



(Ao-Ai)^ 

n 



1/2 



me as 

/Is, (73) 



where tmeas is the time of measurement. This time should 
clearly be much shorter than the relaxation time r^^^ 
(evaluated as « 0.1s, see Table 1) and exceed the rise 
time of the oscillations in the tank trisc ~ 1 /^s (the lat- 
ter condition nicely agrees with the requirement SNR 
> 1). Thus, a drive pulse duration of ^ 10 /iS ensur- 
ing imcas ~ 10 /is seems to be a good choice as it yields 
the sufficiently high value of SNRw 6. The latter (as well 
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resonance tank circuit introduces only small noise at the 
high transition frequency of the qubit, fl^ luq, the rate 
of relaxation can also be made sufficiently small. Ap- 
plying an rf drive pulse of limited span allows a readout 
of the qubit state in the regimes of single and repeated 
measurements. 

Other problems in engineering Josephson qubits with 
weak decoherence are the 1// noise of critical currents of 
Josephson junctions js^ and the 1/f background noise 
coupled to the charge variable [s^ which have not been 
addressed in this paper but are equally important. Hope- 
fully, in the given system these effects might not be as 
strong as in "traditional" tunnel-junction devices like 
dc SQUIDs and single electron transistors operating at 
nonzero voltage bias. Due to perfect decoupling of the 
superconducting loop with the single charge transistor 
from dc bias lines and due to the coherent nature of 
the Josephson current in the zero voltage bias regime, 
one could expect a minor backaction effect of the zero- 
bias operating transistor on its critical current noise and 
charge noise which dramatically depend on the current 
fed (see, for example, Ref. [STfl. 



FIG. 6: The principle of narrow-band radio-frequency read- 
out of the qubit. (a) The resonance curves of the uncoupled 
(dotted line) and coupled to the qubit tank circuit biased at 
operation point A in the excited state (dashed line) and in 
the ground state (solid line), (b) Driving pulse applied to the 
tank circuit (top curve) and the response signal of the tank 
in resonance (the ground qubit state, bottom curve) and off- 
resonance (excited state, middle curve). A smooth envelope 
of the driving pulse is used to suppress transient oscillations 
and has a small effect on the rise time of the response signal. 
For clarity the curves are shifted vertically. 
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as the quantum quality factor Qip) can be substantially 
improved by using a SQUlD-based low-noise amplifier 



CONCLUSION 

We have demonstrated that both manipulation and 
readout of the charge-phase qubit coupled to a tank cir- 
cuit is, in principle, possible. More specifically, the deco- 
herence effect of the electric and magnetic control lines 
can seemingly be minimized by reducing coupling to the 
qubit. The readout system based on the parallel tank 
circuit and cold amplifier can ensure sufficiently weak 
dephasing in the regime without rf drive. The dephas- 
ing rate strongly depends on the accuracy of adjusting 
the offset flux bias <i>m = or = $o/2 correspond- 
ing to operation in the magic points. High symmetry of 
the Josephson junction parameters may further improve 
the coherence characteristics of the qubit. Since the LC 
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